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1. Introduction 

The Y-system for the full spectrum of energies/dimensions in the planar AdSs/CFT4 sys- 
tem conjectured in M has passed a few important tests. It was re-derived and better 
understood within the TBA approach [Q, ||, [|] and successfully tested in the weak cou- 
pling by comparison with the perturbative expansion in N=4 SYM theory up to 4-5 loops 
H ^, p], pi M. Remarkably, the very same Y-system was shown to be responsible for the 
spectrum of /3-deformed N = 1 case in [l(J where the perturbative results || |ll]] were 
reproduced up to 11 loops . At strong coupling the Y-system perfectly reproduces the re- 
sults of quasi-classical quantization of highly excited states of the superstring in 5^ sector, 
in the regime where there is no way to ignore the finite size effects [12], demonstrating 
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Konishi state 




Figure 1: Numerical results of solution of the AdS/CFT Y-system (in integral, TBA form) for 
Konishi dimension as a function of 't Hooft coupling A (from Jl5| and fl6|). 

that the formidable wrapping problem finds its successful resolution within the Y-system. 
Furthermore, some of these results were extended to the generic finite-gap string state in 
|1|] and to the ABJM model in 0. 

The numerical calculation of the Konishi dimension from the Y-system, combined with 



the TBA approach in [15], has provided the data covering a range of values of the 't Hooft 



coupling enough to confirm the leading strong coupling asymptotics 2A 1//4 obtained on 



the string side of the duality [17] and to predict, with a reasonable accuracy, the next, 
subleading correction as being 1.99 • A -1 / 4 . 1 A more extensive numerical study recently 
done in ]l6[ allowed to extend the range of the interval of 't Hooft couplings in the strong 
coupling regime by more than 3 times. This new data from [|l(|] perfectly confirms old 
results: the new points follow the fitting curve from fllBH , within the accuracy margins 2 
(see Fig.g). 

By now it is fair to say that the Y-system is the correct framework to study the 
spectrum of this important duality. However there are still many open problems. The main 
problem which has been preventing us from applying the Y-system to more complicated 
states is to convert it into a system of integral equations suitable for numerical/analytical 
study at intermediate couplings. This problem was only solved for some simple states like 
the Konishi state. It is important to understand the properties of a general solution of 



1 It is tempting to think that it is exactly 2 • A -1 / 4 . We hope that this coefficient will be eventually 
compared with the direct worldsheet 2-loop calculations. 

2 This means that the fitting function presented in |l5[ describes all available data concerning Konishi 
anomalous dimension, including Jl6| , with the accuracy ~ 0.001, which was the target accuracy for At 
the same time, the claimed accuracy of the recent numerics is higher. 
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Y-system describing the physical states of AdS/CFT correspondence. 

The AdS/CFT Y-system is an infinite set of functional equations on the functions of 
a spectral parameter, related to the Hirota bi-linear difference equation (T-system) well 
known in quantum integrability. It differs from the Y-systems for the previously known 
quantum integrable models (ID spin chains, sigma-models etc.) by the specific, so called 
T-hook boundary conditions in the representation space (w.r.t. the discrete indices living 
on a 2D lattice (a, s) G Z x Z) and complicated analyticity properties w.r.t. the spectral 
parameter. 

On the other hand, it is known from our experience with some relativistic sigma models 
and spin chains that the same quantum systems described at a finite volume by a Y- 
system may obey a rather different, finite set of non- linear integral equations (which we 
will abbreviate here as FiNLIE). The first example of such equations, called DdV equations, 
were given by Destri and deVega for the Thirring model [18] followed by a considerable 
activity in this direction. A new approach for the search of such FiNLIE for the integrable 
sigma-models based on the integrable Hirota dynamics of the Y-system, together with a 
few simple assumptions on its analyticity structure, was proposed in [19] and developed in 
|20| ]. It relies on the fact that the Y-system, or the underlying T-system, can be solved in 
terms of a finite number of Q-functions, analogs of those introduced by Baxter for the XXZ 
chain. The T-functions can be represented in terms of finite determinants (Wronskians 3 ) 
of those Q-functions and, knowing the analyticity properties of Q's, one can then write a 
FiNLIE solving the finite size problem. 

Such a FiNLIE would be a very welcome progress for the study of the spectrum of 
AdSs/CFT4 system. It would provide us not only with a more efficient analytic and nu- 
merical tool for the study of this complex model but also most probably give some insight 
into the structure of this duality. However, due to the complexity of the T-hook boundary 
conditions and of the analyticity properties of Y-functions, this FiNLIE remains unknown. 
We propose in this paper an important, in our opinion, step towards the derivation of such 
equations by writing a general Wronskian-type solution of the corresponding T-system for 
the T-hook boundary conditions. This solution will be a natural generalization of our ex- 
plicit general solution, found by three of the current authors [|jj| for a simplified T-system 4 , 
relevant to the quasiclassical limit of the string sigma-model, in terms of characters of spe- 
cific infinite-dimensional representations of U(2, 2|4), parameterized by 8 eigenvalues of 
an arbitrary group element. We will also demonstrate this construction in the asymp- 
totic, large size limit (of long SYM operators) when the asymptotic Bethe ansatz (ABA) 
f22| . |23[| is applicable. We will also discuss the analyticity and symmetry properties of these 
Q-functions. The relation of the character solution to the quasi-classical limit of the su- 
perstring on the AdS§ x S 5 background discussed in [13] will serve as an important source 
of inspiration for that. 

The T-system first appeared in the context of the quantum integrable systems for the 
gl(2) algebra in [24], and generalized to gl(N) in [25|. Wronskian determinant solutions 



3 To be precise, they are variations of a discrete analogue of Wronskian called "Casoratian" . 
4 called in mathematical literature the Q-system, with no explicit spectral parameter dependence (see 
for examples, pl| ] and references therein) 
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of T-system were introduced in for N = 2 and in [27] for any N, where the finite 
dimensional representations of gl(N) symmetry impose the boundary conditions in the 
semi-infinite strip of the width N in the (a, ^-representation space. The Wronskian solu- 
tions were generalized to the supersymmetric gl(M\N) algebras for the T-system 5 with the 
(M\N) fat hook boundary conditions for (M,N) = (2, 1) in g |l| and for any (M,N) 
in p2]] . On the way to constructing similar solutions for the T-hook boundary conditions 
corresponding to the superconformal PSU{2,2\4) symmetry of the model we will first re- 
mind the "classical" solution for the U(2, 2|4) characters [l^, [l3|] and then discuss the most 
general form of the so called TQ-relations relating T- and Q-functions (analogue of Bax- 
ter's famous relation for the XXZ model), in the form of the so called generating functional 
|1, for finite M, |9, 



], and even for the infinite-dimensional [34] representations of 
gl(M\N), as well as the QQ-relations among the Q-functions, especially efficient for the 
super-algebras [jjj (see also earlier papers |56|, The Backlund solution of Hirota 

equation for the gl(M) algebras [27] and gl(M\N) super-algebras, for the fat hook bound- 
ary conditions in the representational (a, s) space |33|], was generalized to the case of a 
general T-hook in [37]. The latter can be conveniently rewritten through the generating 
functional 13] further used in this paper. 

The Wronskian solution for the AdS/CFT Y-system which we are proposing in this 
paper summarizes in the most explicit and concise way all these developments. Its main 



advantage w.r.t. the previous solutions |37, 13] of the Y-system (and the associated T- 
system) in the T-hook is its absolutely explicit form in terms of the Wronskian determinants 
of a finite number of Baxter's Q-functions (7 independent Q-functions) which does not 
include any infinite sums or integral operators. We consider this Wronskian representation 
as a good starting point for trying to derive a FiNLIE system describing the spectrum of 
the planar AdSs/CFT,^ 



2. F-system and T-system for the spectrum of AdS/CFT 

In this section we will remind the Y-system for the spectrum of AdSs/CFT,^ point out its 
symmetries and discuss its analytic properties. 



a 



O O 



s 




O O 



Figure 2: T-shaped "fat hook" (T-hook) uniting two SU(2|2) fat hooks, see J| for this T-hook 
and its generalization fj7|. 

5 T-system for this case appeared in feq, E9|. 
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The problem of the spectrum of the closed superstring sigma model on the AdSs x S* 5 
background, similarly to all known integrable sigma models on a space-time cylinder with 
the global gl(M\N) type symmetry, can be reduced to the Y-system || 

y+y- (l+yg, 8 +l)(l+yq, 8 -l) , 21 v 

«,. a,s (1 + i/ Fa+ls ) (1 + 1/F a _ 1>s ) 1 • > 

where Y a ^ s are functions of the spectral parameter u, defined on the visible nodes of the two- 
dimensional integer (a, s)-lattice presented on the Fig.Q(right). By the subscripts of any 
function f(u) of u, we denote the shifts of the spectral parameter by ±|: = /(u± |). 6 

What distinguishes the Y-systems for various sigma models is: 1) the boundary con- 
ditions on the a, s lattice which are mostly defined by the symmetry of the model 2) the 
analytic properties w.r.t. the spectral parameter u, only partially constrained by the sym- 
metry and greatly depending on the physical properties of the model. 

The Y-system for AdS/CFT, as any Y-system, is directly related to the so called 
T-system. Namely, defining 

v _ ?q,s+l^a,s-l ,„ „n 

X a,S — rp rp 

we rewrite ( |2.l| ) in terms of the bilinear finite difference Hirota equation 

Ta, s T a , s = Io+i )S T _i )S + T aiS+ iT ajS _i . (2.3) 
For our string sigma model the boundary conditions on the a, s lattice are constrained 



by the superconformal PSU(2, 2|4) symmetry [13[ and impose that the T-functions, are 
nonzero only inside the so called T-hook i.e. only on the nodes of the visible part of 
the (a, s)-lattice of the Fig.^(left). The rest of the T-functions are zero: 

T a>s = for a < OU (a > 2n |s| > 2) . (2.4) 

These boundary conditions for the functions T as agree with the above mentioned 
boundary conditions for the functions Y a>s of the Fig.|2|(right) if we take Y-functions to 
be zero on the vertical boundaries, and to be infinite on the horizontal boundaries on the 
Fig.^(left). Note that this leaves an ambiguity about the values at the corner nodes Y2,±2- 
A more careful analysis shows that equations for Yi,±2 are more complicated and cannot 
be written in a "local" functional form. The two missing equations ( |2.11| ) do not have the 
standard Y-system form and can be borrowed from the TBA approach [Q]. However, we 
believe that all the functions T ajS satisfy the standard Hirota equation ( p.3[ ) Q. In that 
sense, the T-system looks more universal, and in many cases more convenient than the 
Y-system. 

Notice that the parameterization in terms of T a)S is not unique. The T-system is 
invariant under the gauge transformation: 

T a ,s -> 9 [ r +s] 9 l r S] 9t a] 9[- a - S] T a ,s (2.5) 



6 and in general / [±s] = f(u ± f ). 
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and thus the T- functions are defined up to 4 arbitrary gauge functions g a (u). The physical 
quantities are computed in terms of the gauge invariant Y-functions. 

The analytic structure of the AdS / CFT Y-system is inherited to a great extent from 
the dispersion relations for the elementary physical excitations - the magnons on the in- 
finite spin chain representing a SYM operator, or its AdS dual on the string side. The 
dispersion relation between the energy and the momentum of such solitary excitations [38] 
is conveniently parameterized 

is related to the 't Hooft's coupling A 



9| by the so called Zhukovsky map ^ = x + 1/x, where g 



9ymNI 



*> 9 = & 



p(u) 



TlOg 

% 



x + (u) 
x~(u) 



e(u) = 1 + 



2ig 



2ig 



~(u) x (u) 



(2.6) 



The inverse map is double valued and we have to distinguish two branches - physical and 
mirror (related to the exchange of time and space coordinates on the world sheet cylinder 
li): 



u) 



9 



+ 



+ 2 



9 



9 



+ i\ 



.9 Z 



(2.7) 



In the physical branch, the finite cut, by definition, connects two branch points u = ±2g 
whether as in the mirror branch the cut connecting them passes through the infinity. 

The u-parameterization is distinguished by the fact that the fusion of elementary exci- 
tations into various bound states is especially simple in the complex u-plane: the rapidities 
of the constituents of a bound state are spaced by the integers of i. For example, the bound 
states for the energy carrying magnons mentioned above have the energy and momentum 



m 



expp a (u) =J„o expp(u), exp e a {u) = T a o exp e(u) 
where we defined the fusion operator 



?a f( u ) = < 



[qj — 1 



11 f(u + ij) , a>0 
1 , a = 



a I — 1 

*=-H- 



H-i 



11 !//(« + ij) , a<0 



lal — 1 



This gives for the dispersion relation of the bound states 



(2- 



(2- 



1 / x^\ 2ig 2ig 

Pa{u) = T log Uh ' £a(M) = a + ^ ~ 



(2.10) 



The Y-functions as functions of u should inherit the multi-valuedness of the map (^jj) • 
Most of our experience on their analyticity properties comes from the asymptotic Bethe 



ansatz (ABA) [p2| , 23] corresponding to the limit of very long operators L —¥ oo and from 
TB A equations for the excited states [0, 0, 0] . The ABA limit of Y-system found in Q 
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shows that the Y- functions have branch points at u = ±2g + ^ for various n's (re G Z). 
Similarly to the above definition of x mir we define the "mirror" sheet of the Y-functions 
with the cuts going through infinity, parallel to the real axis. We can find from the study 
of the ABA limit (see the section ||) that Yi jS | s>2 has 4 branch cuts at Im u = ^ 1 , while 
^a,±i| a >2 nas ^ branch cuts at Im u = 2 1 an d ^,±1, ^2,±2 have 3 cuts at Im u = 0, ±1. 
At finite size other cuts can appear, but we expect the following analyticity conditions to 
hold anyway p3[ |: 

Isl — 1 \s\ — 1 

1. Yi t ± s have no branch cuts inside the strip — ■ u 2 — < Im u < U 2~; 

2. Y a: ±i have no branch cuts inside the strip — <Imu< 

3. Y a fl have no branch cuts inside the strip — § < Im u < |; 
4- ii,±l) ^2,±2 have a cut on the real axes such that 

Y l!±1 (u + iQ)Y 2 , ±2 (u-iO) = 1, u G (-00, -2g] U [2 5 , 00); (2.11) 



5. y ajS obtained from the generating functional should be real functions in "mirror" 
kinematics. 

The Y-system should be satisfied for the mirror branch of Y a ^ s described above (i.e. 
for the cuts chosen to go through infinity). We also define the "physical" branch of Y\q 
as an analytic continuation through the first cut above the real axis and then back to the 
real axis (as proposed in |i~5[| ), and denote the continued function as Y^q. These properties 
are very similar to the ones used in [19| to convert the su(2) principal chiral filed theory 
Y-system into the corresponding FiNLIE. 

Once the appropriate solution of the Y-system for a given physical state of the 
AdSs/CFT4 system is found its energy is given by the following formula 

£ = E e i h Ki) + E / S %r lo § (! + Y M) , (2.12) 

j a=l 

where we also have to impose the Bethe ansatz equation for the Bethe roots (the rapidities 
of physical excitations) [§, |l5[ 

n P oK,) = -l. (2.13) 

the exact Bethe equations for the auxiliary roots should come from the condition of pole 
cancellation. 

In the subsequent sections we first recall the solution of the Y-system in the classical 
large A limit given in terms of the characters and then describe our construction for the 
general quantum solution. 
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3. SU(2, 2|4) character solution of the Y-system and classical limit 

In the classical limit when A is large the shifts by ±i/2 in the (|2.3|) become irrelevant (see 



12] for more details) and the functional T-system reduces to an algebraic set of equations 



called in the mathematical literature the Q-system 7 : 

(T a ,s) 2 = Io+i )S r a _x lS + T a)S+ iT a>s _i . (3.1) 

In the paper [l^] three of the current authors clarified the group theoretical meaning 
of the AdS/CFT y-system, its relation to the characters of irreducible representations 
of the SU(2,2\4) symmetry, related to the superconformal PSU(2,2\4) symmetry of the 
model, and their relation to the classical limit, extending some of the results of Jl^] to all 
sectors. We will briefly remind in this subsection the basic results concerning the explicit 
construction of the SU(2, 2|4) characters for the unitary representations, in terms of finite 
determinants, in full similarity to the 1-st Weyl formula known for the compact represen- 
tations of GL{N). In the next chapter we will show how to generalize these formulas to 
the quantum solution of the T-system ( |2.3| ) , and to find the general explicit solution of the 
underlying y-system and T-system in terms of finite determinants, called Wronskians, pa- 
rameterized by a finite number of Baxter's Q-functions. We hope to apply these results in 
the future for the construction of a finite system of non-linear integral equations (FiNLIE) 
describing the full spectrum of the planar AdSs/CFT^ 

3.1 The character solution of the simplified AdS/CFT Y-system in T-hook 

The generating function of U(2, 2 1 4) characters of "symmetric" representations can be 
represented as 

(1-^(1-^) (l-y^jl-yjt) 

W ^ k) - (1-^(1-^) X (1-X 3 t)(l-X4t) (3 ' 2) 

where (xi, ...,x^\y^... , j/|) are the eigenvalues of a group element h € U(2, 2|4). Here 
and below we use the hats over indexes to indicate their "fermionic" grading, as opposed 
to the "bosonic" grading of the rest of indices. 

The first and the second factors in the r.h.s. may be attributed to the right and left 
U(2\2)r^l subgroups. The characters of "symmetric" representations are generated by 

dt WM4(t; h) 

JC 27TZ 

where the integration contour C encircle t = together with the poles — , — correspond- 
ing to the subgroup U(2\2)l, leaving outside the poles corresponding to the first 

subgroup U (2\2)r. Note that s can here be positive as well as negative: — oo < s < oo and 
the corresponding irreps, first constructed in [44, 45], are infinite-dimensional: T-^'^f/i] are 



Tif\h}=$ — rr"V -i, (3.3) 



not polynomials of anymore, unlike the compact representations of U(M\N). The 

rest of the T Q)S (/i) can be restored by means of the Jacobi-Trudi type formula: 

T a s = det Ti s+i-j ■ (3.4) 



7 It should not be confused with the Baxter's Q-functions considered below. 
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It is easy to see (at least on Mathematica, using the code given in jDJ) that T a ^ s 7^ only 
for (a, s) G T where by T we denote the T-hook drawn on FigjJ(left). In [13] the explicit 
expressions for all these characters in terms of 2 x 2 and 4x4 determinants was found: 



det ( ^+y- 4 - (a+2)e -+ 2 



( — \ \(a+l)s / £324 

K'li^mJ dct(s^' 0+ y~ 4 ~ (0+2)8j < 0+2 ) 



a VJ 1 -^-) a-j+(«-a)(i-e io ) 

dot Z; £• 



de t (Z (1 - £ '- 0) a; 2 - +(0 - 2)(1 - fl - 0, 



(3.5) 



where 5? 



fa;-^3)fa-^4) ^ 
fe-a:i)fe-^2) ' * 



The 



I 1 , 3 > s 
{ , j<s 

other T's can be obtained using the wing-exchange symmetry which is related to an outer 
automorphism of the Dynkin diagram of g[(4|4) 



T a}S (xx,...,x4\yi,...,yi) 



1 

X4 



T a —s I ' • • • ' 



1 

Xi 



1 



\x1X2X3Xi J 

Another important property, under the rescaling of eigenvalues, reads 

T a , s (axi,ayj) = a as T a)S (xi,yj) 



1 



(3.6) 



(3.7) 



which implies that the gauge invariant quantities are invariant under the P-transformation 
from PSU(2, 2|4). Note that the Weyl group symmetry w.r.t. permutations of (x±, X2, £3, £4) 
is broken: these characters are symmetric only w.r.t. x\ O X2 and X3 o X4. 

As we mentioned in the beginning of the section this character solution in the T-hook 
is tightly related to the classical limit of the AdSs x CFT4 Y-system. In this limit the string 
is long, the AdS/CFT coupling is large and the natural scale of the spectral parameter is 



u ~ g. We see that we can neglect the shifts w.r.t. the spectral parameter in (2T) and 
(only a slow parametric dependence of the group element h(u) is left) and Hirota equations 



takes a simplified form, called Q-system (3.1). Note that ( |3.5|) represents the most general 
solution of this equation with the T-hook boundary conditions and is parameterized by 8 
eigenvalues 8 . 

It was shown in |l3| that the characters ( |3.5| ) give the classical limit of the AdS/CFT 
T-system (or Y-system) if the group element h is simply the monodromy matrix f2(ii) of 
the classical finite gap solution of Metsaev-Tseytlin superstring [46, 48, 49. [5(| 



T a , s = Str aiS ft . (3.8) 

Thus the solution of Hirota equation can be expressed solely in terms of the eigenvalues 
of the classical monodromy matrix. The eigenvalues as functions of the spectral parameter 
represent 8 sheets of the finite gap algebraic curve and the generating function fl3.2|) can 
be viewed as its spectral super-determinant. 

8 One eigenvalue can be always rescaled to unity leaving us with only 7 gauge independent parameters. 
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3.2 Quantization of the classical algebraic curve 

For the states with large length L the spectral equations should simplify to ABA equations. 
These states are described within the ABA approach by the configurations of the Bethe 
roots. A map between the configurations of the Bethe roots and the classical algebraic 
curve was understood in p7], 51, 34, ^3, [52|, 53]. The basic idea is that the branch cuts of 



the classical curve are the cuts on the complex plane where the Bethe roots are densely 
distributed. Technically one can write the classical quasi-momenta in terms of the densities 
of the Bethe roots or equivalently resolvents: 

Ka 2 1 

H * = Y, -T-T ' S ^ x ) = H *Q-/ X ) ■ (3-9) 

where a = 1, . . . , 7 is a type of the Bethe root. The eigenvalues of the monodromy matrix 
are related to the quasi-momenta by the formulas (see the notations in [EH, 54]) 



(3.10) 



Vi = e -<fc = exp - xH x - iS 3 + iH A ] 

Xl = e -<Pi = exp Lj WtWQi - iH x + iH 2 + iH 2 - iH 3 
x 2 = e -*P2 = exp (-i Lx/( ^ ) _ + l lQl - iH 2 + iH 3 + m - iH 2 
y - = e -& = exp + iHs - iH A + iH x 

y . = e -m = exp _ m + m _ iSr < 

X3 = e-ws = exp Ui Lx ^)_+iQi + iH& _ _ iSj + iH 6 

X4 = = exp (+i W gy gl + iH 7 - iH 6 - iH & + iH 5 
Vl = e -ip* = exp ( +i W(2g)+ig a » + iHr + iBb _ iS 

Where Qi = ^ - ^ , Q2 = Yli 7T - Notice that only two subsequent sheets 

share the same resolvent. This discretization is not unique. In particular, a configuration 
of Bethe roots could be mapped to a different one using the so-called fermionic and bosonic 
dualities, which correspond to a reshuffling of the sheets of the Riemann surface (e.g., see 
fl~3[| ). One can make a permutation of sheets and discretize the quasi-momenta in terms 
of a new set of Bethe roots defining new resolvents (see Fig.|3|). The cuts connecting some 
non-neighboring sheets can cross on the way the other sheets due to a mechanism of "stuck 
formation" of various types of Bethe roots |54| (strictly speaking the stacks are only formed 
for small enough filling fractions and nonzero twists |55|] ). 

We will see that it is useful not to restrict ourselves to a single possible discretization 
(related to the quasi-classical quantization) of the curve but rather to consider all possible 
permutations of the sheets. These permutations can be denoted by sequences of indices 
labeling the sheets. In these notations, we denote the chosen ordering (grading) of sheets 
as I 1223344. 

The grading 11223344 has some obvious advantages in the ABA limit and we will use 
it as a canonical one. 
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Figure 3: One can quantize the classical curve by choosing different orders of the sheets. Once 
an order (called path) is fixed the Q functions will correspond to the cuts crossing the line between 
two neighbor sheets. A Q function is not affected by the permutations of sheets above and below 
the line. Thus the Q function is uniquely determined by the subset of the sheets above the line and 
we denote it as Qi 12 23- F° r more formal definitions see section 4.3. 



The main motivation of this paper is the generalization of this solution to the full 
quantum case of the AdSs/CFT^ Y-system valid for any operator and any A. In the next 
section we show that at least the first part of this program - the general solution of the 
T-system (|2.3| ) in the T-hook, in terms of a finite number of functions of u — the Baxter's 
Q-functions - can be fulfilled in a rather explicit form. The next step - finding the general 
analyticity properties of these functions and constructing the corresponding FiNLIE, will 
be the subject of a future work. 

4. Solution of the full Hirota equation in T-hook 

Our goal here will be the construction of a general solution of Hirota equation ( |2.3| ) for the 
(2|4|2) T-hook, first in terms of the generating functional (GF) and then in explicit and 
finite form, in terms of the Wronskians - the finite determinants of Q-functions. 

4.1 Generating functional for the T-hook 

For given T\ s (u), —oo < s < oo, the general solution of (|2.3| ) is given by the Bazhanov- 
Reshetikhin (BR) type formula 9 

T - det rA-a+k+j-l] (4 1) 

As was shown in |28|, p9L |33[| , to solve Hirota equation in a (M\N) fat hook the "sym- 
metric" T-functions Ti iS (u) should be generated by a special generating functional. In 



9 The original BR formula was written for the fusion in spin chains. It is a generalization of the Jacobi- 
Trudi formula (3.4) for the characters. Here we view as a general solution of Hirota equation on the 



upper part of the (a, s)-lattice, with the Dirichlet boundary conditions: Tq ]S — 1 and Ti jS - fixed. 
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particular, for the (2(2) fat hook, in the grading (1|2|1) corresponding to the Kac-Dynkin 
diagram ®— 0— <S>, or to the {yj|cci, a^lygi or dering of the eigenvalues of a twist parameter 
g G GL(2|2), we have the following generating functional 



00 

W = U t Uri Ur^ U ix ^ = £ D s T hs D s (4.2) 

S = — OO 

where 

U A = (1 - D XA D), (D = e-^) (4.3) 

with A C Io = {1|1,2|2} being an arbitrary subset of the full set Iq, while 
ix. i( u )IXi i( u ),Xii 2 ( u ) IXii2 2( M )i are arbitrary functions of the spectral parameter u, re- 
placing the group element eigenvalues X2I2/2} °^ the character generating function 
( |3.2| ), whereas the generating parameter t is replaced by the shift operator D = e~z du . As 
was proposed in |J2| in the case of a more general superalgebra 10 gl(M\N) (and demon- 
strated on the Fig.^), the labeling here follows the pattern of the Dynkin diagram in the 
grading <^-&- £g>: we start from the l.h.s. of this diagram with A = (, 1) for the first factor 
in the l.h.s. of ([4.2;); then for the label A of the second factor, one moves the comma one 
step to the right and adds new index after the comma, making A = (l,i). The choice 
of the index is related to the grading: when moving to the right of the Dynkin diagram 
we cross the fermionic node which means that we chose to add an entry from / with a 
different grading than the previous i, say 1, which gives for the second factor A = (1,1); 
then we cross the bosonic node, which means that we should add an entree with the same 
grading as the last one, i.e. 2, which gives for the label of the third factor A = (11, 2), and 
finally, crossing the last fermionic node we get for the label of the last factor A = (112, 2). 
This means that the subset A represents a "path" by which it was reached starting from 
the l.h.s. of the Dynkin diagram. This will be the general rule of for more complicated 
algebras, and in particular u(2,2|4). 11 

To calculate T a ^ s (u) in terms of these functions we formally expand the l.h.s. in powers 
of the shift operator D = e~2° u and compare the coefficients with the r.h.s. The formal 



proof can be found in [28, 29, but it is easy to convince oneself on Mathematica that 
the T-functions generated from Ti jS by ( |4.1[ ) are zero outside of the (2|2) fat hook. 

What is the solution for the T-hook of u(2, 2|4)? We can read it off from the generating 



function for the characters of it(2,2|4) given by eqs.( |3.2|) -(3.3). Note that the integration 
contour prescription in ( |3.3| ) means that we can generate the u(2,2|4) characters by ex- 
panding the first factor in the r.h.s. of (3.2) in powers of t, and the second factor - in 



powers of t, and then extract from the product Ti jS as a coefficient in front of the power 
t s . The result is given in terms of infinite sums, and not polynomials in the eigenvalues, 



signalling that we deal with infinite-dimensional irreps of u(2,2\4) p4|, p4|, 45]. 



10 It is important to note that the Baxter Q-operators are essentially independent of the shape of the 
"hook" . They are fixed by a certain oscillator algebra p6| , |56| [57| |3l], |H^] , at least for the models with the 
Yangian or the quantum affine algebra symmetries. This allows to apply the formalism for the (M|7V)-hook 
directly to the construction of the Wronskian solution for the T-hook in question. 

11 We will see later that not all the paths give inequivalent \A,i functions. Namely, = Xi',i' iff I an d 
I' differ only by the permutation of their indices, while i — i' . 
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Similarly to the case gl{2\2) described above, we can try as a quantum generaliza- 
tion of equation ( [2 .3; ) for the (2|4|2) T-hook, for the particular grading of the eigenvalues 
{yi\xi, X2\y2, 2/3 1 , £4(2/4} corresponding to the Kac-Dynkin diagram 0— 0— 0— 0— 0— 0— 0, 
the following generating functional p4"[ 



W = [V X Vf 1 V 2 l Vi] + x [V$ Vf 1 V4 1 Vi] = J2 D'T Xt ,D' (4.4) 

S = — OO 

where we introduced, to make the formula less bulky, the notations Vj = 1 — D£jD with 

fi = Xi & = 6 = Xii, 2 £2 = Xii 2) 2 

= ^1122,3 ^3 = Xil223,3 ^ 4 = ^il2233,4 ^4 = ^il22334,4 • (4«5) 

We label the functions £j only by the last index in the subset A of xa having in mind 
the particular grading, or nesting path. By definition, we expand inside [. . . ] + and [...]_ 
w.r.t. the positive and negative powers of D respectively and then calculate the functions 



Xi jS , — oo < s < oo, comparing the powers of D on both sides of the eq.(4.4). To restore 
the rest of T a s we can use again (|4.1|) . 



The formal proof 12 of the fact that (4.1), (4. 4) represents the complete solution of Hirota 



equation within the T-hook will be published elsewhere, but, again, it is easy to convince 
oneself in the correctness of this formula on Mathematica. 



The representation (4.4) is already a considerable advance w.r.t. the original Y-system 
(2.1) since its solution is now parameterized in terms of a finite number of functions £j. A 
solution of Hirota equation is now parameterized by 8 arbitrary independent functions 
entering the generating functional (<L4) 13 . Moreover in this form the passage to the general 



strong coupling solution found in [13] is straightforward - in the strong coupling limit the 
shift operator D can be replaced by a formal scalar expansion parameter t whereas should 
become the eigenvalues of the classical monodromy matrix. This solution is quite useful 
for various applications. Note however that as the result we get infinite series, signaling 
that we deal with infinite dimensional representations of U(2,2\4). In the next section, we 
will show that these infinite series can be "miraculously" converted to the explicit finite 
dimensional Wronskian determinants. 

4.2 QQ-relations 

To find the Wronskian solution we have to choose a good parameterization for all these 



8 functions (|4.5[). The formalism of p8l 29, 34, 33, p2| aimed at the derivation of Bethe 



equations as a condition of analyticity, i.e. polynomiality, of the T-functions, as well as 
the form of the classical finite gap solution ( 3.10|) suggests that the best parameterization 



would be in terms of the Q-functions (a-la Baxter). As it was explained in sections 3^ and 
4.1 and on the Fig.0, all the 2 s Q-functions can be labeled by all possible subsets from the 



full set I = {1,2, 3,4(1,2,3,4}. 



12 The Backlund procedure of [ p3[ |37| can be used for it. 
ls One of them can be removed by a gauge. 
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Figure 4: Hasse diagram for gl(2\2) (cf. |32|): To construct the generating functional (4.2) we 
move along any path starting at the node and ending at the node 1212. Each line corresponds to 



adding the next [/-factor (or V-factor in (4.4) for the (2, 2|4) case) with one extra index. The dashed 
blue (resp. dotted red) lines stand for adding a "bosonic" (resp. fermionic) index, and each local 
change of the nesting path at some rectangular facet of the diagram gives rise to the QQ-relations 
( |4.7| ), (p~8|). The equation (hS) (resp. (17)) corresponds to the facets having 2 "bosonic" lines and 
2 "fermionic" lines (resp to the facets having all 4 "bosonic" or all 4 "fermionic" lines). 



The monomials of the generating functional can be conveniently parameterized in terms 
of 8 Baxter's Q-functions as follows: 



Q Qf ' 

Q7~.. -Q+ 

11223344 1122334 
Qll223344Qi 12 2334 



Cl 
£4 



Q7Q+ + 
^1 ^11 

QT Q7".. 
1122334 112233 

Q il22334 Q 112233 



6 
£3 



Q7-Q+ 
11 112 

^J- 1 112 



J.J.22.',:- 1 , 11223 
12233^ii223 



£2 
£ 3 



Q+ Q7". 
112 1122 

Q ri2 Q 1122 ' 

Q7 .Q++. 
11223 1122 

Q il223 Q il22 



(4.6) 



The way we perform this parameterization for the functions xa (see the definition ( |4.5| )) 
in terms of the Q-functions is completely denned by the chosen nesting (grading) path and 



is simply given by xi 



[k] 



last Q in the denominator of xi > a 



15 



where Qj in the denominator is the same as the 
= ±1 is chosen so that the first ratio in xi,j is the 



same as the last ratio in xi when we change the grading passing from (I, i) to and 
a is opposite otherwise. All such nesting paths for the case of su{2\2) are shown on the 
Hasse diagram in Fig.|j. 

All that means that the labeling of Q-functions follows the nesting path. However, it is 
clear from the Fig]3| and the explanations in sec.[T^ that these Q-functions do not depend 
on the order of indices in the path up to a sign: a permutation of indices multiplies the 
Q-function by the signature of the permutation 16 . 

Moreover, all these 2 8 = 256 Q-functions are not independent and can be expressed 
through a chosen basis of 8 Q-functions. The rest of them are related to those 8 by the so 
called QQ-relations ||, p], H H © HI 0: 17 

QiQi,ij = Qfe 1 Q^ 1 -Qfc Pll Q! p ; 1 for p l = Pj , (4.7) 



To understand what is the first and the last factor, the second line in (4.6) should be read from right 
to left, which makes the number of indices increase continuously. 

15 which helps to cancel the poles given by zeros of Q-functions, if there are any, in all T-functions by 
imposing the Bethe ansatz equations on their positions. 

16 The ± sign coming form this signature is irrelevant in the Xa,i functions, where these signs are canceled 
in the ratios of Q functions. The statement of the footnote |ll| is a consequence of that. 

17 Of course, any linear combination of solutions of the Baxter equation is also a solution. Besides that, 



Ql,iQi,j = Qj Q/Ji - Qj Q/7y for Pi = ~Pr (4-8) 



where Pi = +1, i £{1,2,3,4} and Pt = -1, i£ {1,2,3,4} (4.9) 



These QQ-relations can be obtained from the fact that the relabeling of ^-functions 



and Q-functions in (4.4) by changing the nesting path does not change the generating 
functional. This change occurs from the commutation of two consecutive [/-factors in the 
generating functional. Following the nesting path should give the same 



generating functional in 4.4 on page 13 as the path — > 



V7! iy ui '= v u JV ^ ( 4 - 10 ) 



which gives both "bosonic" and "fermionic" QQ-relations ( |4.7D and 

Let us note that the dependence on the general (u- independent) twist matrix 
g = diag{zj|zi, z%\z§, z^\z^, z^zg} S GL(4|4), useful for some applications, such as the 
/3-deformed version of the AdSs/CFT4 duality, can be easily introduced into the T-hook 
solution ([4.4p by a certain simple rescaling of all Q-functions. 18 



4.3 Explicit Wronskian formulae for T-hook 



Now we can formulate the central result of this paper - the determinant formulae for the 
complete solution of the Hirota equation (T-system) fl2.3| ) within the (2, 2|4) T-hook. They 

by the following three formulas 



are written in a gauge where Qg = l 19 , and the whole set of non-zero T-functions is given 



we have a rotational symmetry GL(4) x GL(4) of QQ-relations. In addition, there are discrete symmetries 
related to the outer automorphisms of pZ(4|4) and the gauge symmetry. However, once we fix a basic set of 
8 arbitrary Q-functions, we will fix the 2 4+4 Q-functions unambiguously, through the QQ-relations. 

18 One can recover the twist {zt} after the following formal replacement: Q^ ' — > aif^Qf\ tL°' s — > 

«1234i234/L7i234 T ™ > ^ = IWi" 1 > /J* 1 = fP - = W^I^^jY^ (<* eqS. (3.26)- 

(3.34) in [B2[). To take a limit Zi —¥ 1 is a non-trivial task, where one have to take into account the 
rotational symmetry of the QQ-relations (see recent papers JgoJ ) . This produces an extra factor of Zj in 
the r.h.s. of each of the eight expressions for £j in ( J4.6| ) . 

19 We will now focus on gauges where Q0 = 1, in this whole article, since Q0 7^ 1 case can be easily derived 
from the formulae of this paper by a gauge transformation. 
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Ti,s 

T 2,s 
Ta,+2 
Ta,+1 

T a ,0 
Ta,-1 
T a ,-2 

T i,- 



+Qf Qp - Q^Q^l ,0+1 



<2 - 2 

s ^ +2 



(4.11) 



V I V m 121 122^122 ^123^123 124^124 J ' 

+q[°Lq [ ~?I - o [a Lo [ - a ! + Q [a Lo [ -?l +Q [a Lo [ - a l - o [a Lo [ - a l + o [a Lo^ a l 

^1212^4343 ^1213 4342 ^1214^4332 ^1223 4341 ^1224^4331 ^1234^4321 

V 7 V 434 434 W 433 433 ^ ^432^432 ^431 431 , 1 

+Q^ sl QS , , < -2 
+Q^ s] Qi sl - Qt^Qf 1 , 



where the bar over the indices denotes the complementary set, for instance 

O •••-21 



<34i234' Note that this solution obeys the following properties: 



(d u - 2id s )T 0tS (u) = 0, 



T 



k,2 



2,fc, 



fc,-2 



T 



2,-fcj 



k > 2. 



(4.12) 



That means that we have imposed three out of four possible gauge constraints ( |2.5[) . 

Using the determinant formulas ( |A.1|) and ( |A.2j ) solving the QQ-relations and the 
Laplace expansion formula for the determinants it is also possible to rewrite T-functions in 
various determinant forms. For example, we have the following 4x4 and 2x2 determinant 



expressions reminding the character formula ([ 
we can represent T a>s as a 4 x 4 determinant 

(_!)«. 



For the upper strip a ^ \s\ of the T-hook 



(^ s »Qil)(^ 

and for the right strip (s ^ a 0) 



Q 



-a\ 
34 



(r\[a.—s-2i+'i}\ 

V^12j Jl^i^2-s,l<j<i 
(r\[— a+s+3-2i]\ 



-1 



° Q 



341234 



Q 



[-s+a-3+2fc]\ 



.7,341234 



l<fc^2-a 



Q 



[s-a+2fe-l] 



lsgjr'SC2 

lsgfc^a 



(4.13) 



(4.14) 



20 The bar denotes the complementary set when I is sorted in increasing order for the canonical ordering 
1<2<3<4<1<2<3<4. For instance, Q^j = Q^i^i- K the set 7 is not ordered this way, one 
has to commute its element first (which introduces a sign), and then to use this definition of the bar. For 
instance = -Q^ = -Q 12 i 234 . 

21 Each term of these formulae should be T-functions related to infinite dimensional representations of 
9 Z(4|4). Let us introduce functions t M = Q^Q^T* 1 - Q^Q^"" 1 and T M = -Q^Q[" s| + Q^Ql" 5 '. Note 
that the above T^J hook coincides with Ti jS for s £ and Ti, s for s 6 which correspond to infinite 

dimensional irreducible representations. Then we find that the T-function related to finite dimensional 
irreducible s-th symmetric tensor representation of pZ(4|4) is a subtraction of these functions: x^~ hook = 
Ti >s — Ti. s |3^ ]. Thus we see that an analytic continuation of Ti jS to positive s should be a T-function for 
an infinite dimensional reducible representations, which contains the finite dimensional irrep. This is similar 
to Bernstein-Gelfand-Gelfand-resolution on the level of the T-functions. 
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and the left strip (— s ^ a ^ 0): 

'„Ls-ra+2ifc-1l\ /„f-.s+a-3+2fcl\ 

(4.15) 



1 a ^121234 



n [s-a+2fc-l]\ / n [-s+a-3+2fc]\ 
S J3^j^ ^I2j,i234 J 3^4 



as 2 x 2 determinants. 

Note that the above T-functions are taken in a different gauge than the T-functions 
of the generating functional. They are related as Taj s = qL s ^Q' ^^^1°], where Q0 = 1 
in our normalization. It is not difficult to see that these formulas solve indeed the Hirota 
equation in T-hook. Indeed, the Wronskian determinants ( 4. 13| ) , ( fl.14] ) and ( 4.15| ) have 



similar structures 22 as the ones solving the Hirota equations in the infinite strips of a sizes 



N = 4,2,2, respectively (see ]27j, and its applications [|2C|]). 

It is easy to see that to saw the three semi-infinite strips together into the full T- 
hook it is enough to prove that the Wronskian solution satisfies four Hirota equations 
corresponding to the nodes (a, s) = (±1, ±2). This can be done straightforwardly. 

There exist many possible determinant representations for the solution. 23 In the next 
subsection we will show, as an example, how to write every Q-function explicitly in terms 
of a basis of 8 Q functions. 

4.4 A basis for Q functions 

The aim of this subsection is to show that it is possible to express explicitly and in a finite 
form all the T-functions of the T-hook in terms of a basis of 8 Q-functions, at least for a 
particular choice of this basis. 

As we already mentioned, all the T-functions can be expressed, as the general solu- 
tion of Hirota bilinear difference equation, through 8 independent functions. Already the 
generating functional (|4.4j ) performs this task but, unfortunately, every T-function is ex- 
pressed only by an infinite series in terms of 8 £- functions ([4.5|), organized in powers of 



wrapping. Another solution involving infinite series was recently proposed in [37] which is 
closely related to the generating functional ( f4.4| ) presented several years ago [34| in a dif- 
ferent context. Their parameterization ( |4,6| ) in terms of 8 particular Q-functions does not 
immediately lead to finite expressions for the T-functions since to express the Q-functions 



of ( 4.11 ) in terms of these 8 Q-functions we need to solve the appropriate QQ-relations 
which, a priori, will lead again to the inversion of some finite difference operators, and 
hence to some infinite products. 

Nevertheless, we will show here that all the Q functions can be explicitly expressed in a 
finite form as Wronskian determinants in terms of a particular basis B\ = {Qi, Q2, Q3, Q4, 
c/j, q^, gg, ^4} where = Qj/Q^. First, the bosonization trick (see Appendix^) says that 



by relabelling all the Q-functions as Qj = Q/ASi defined by (A4), one gets Q-functions 
which obey for any indices (for the gradings with or without hats) only the "bosonic" QQ- 
relation (4/7) with = pj = — 1, Vi,j. In other words, the functions Q don't distinguish 



22 However one has to note that the matrix elements and the prefactors of the Wronskians are quite 
non-trivially related each other by the QQ-relations, in this case. 



23 Some of them will be published elsewhere [ plf . 
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bosonic and fermionic indices, and all the Q-functions with more than one index can be 
expressed through one-index Q-functions by standard bosonic determinants, namely 



Qi! 



=j[-JV-l+2Z] 



(4.16) 



Jjv-i ° Q0 

In this manner, every Q-function can be represented as a determinant of the Q- 
functions of the set B 2 = {Q0, Qi, Q2, Q3, Q4, Qi, Q 2 , Q 3 , Q|}- 

We will show now that we can express all Q-functions through 8 functions of the set 
B\ in an explicit and finite form. On the other hand, £> 2 contains 9 Q-functions, which 
are not all independent. Due to the Pliicker relations among various Q-functions given 
in the Appendix [A], these 9 functions can be expressed as determinants of 8 independent 
functions of the basis B\ as follows 24 : 



Q 
Q b 



Q [j] , 

j=3,l,-l,-3 

W J I lCi«;4; i+b 
j=2,Q,-2 



W / ^0 



l<i<4; i^f 
j=2,0,-2 



1 < b < 4 



(4.17) 



(4.18) 



(4.19) 



where 



Qi = Q 



l<£sg4 
-1,-3,-5,- 



Q-/Q- 



(4.20) 
(Q234, Q134, 



j=- 1,-3,-5,-7 

The relations (p^|) and (|1|), for Q = Q1234 and (Qi, Q 2 , Q 3 , Q 4 ) 
Q1247Q123) are simply the determinant solution of bosonic Pliicker QQ-relations ( |A.1[ ) (in 
the gauge = 1). In the same manner as the bosonization trick, the transformation 1 1— >■ I 
simply exchanges the gradings of all indices 25 : so that the same bosonic determinant gives 
for instance 



det 



Q 



1^3 

J=2,0,-2 



Q 



[j] 



det 



123 



Q±Q: 

i 



and Q 



1^4 
j=3,l,-l,- 



Q 



[.;] 



1234 



Due to the definitions (Qj, Q 2 , Qg, Q^, Q0) 
nants are easily recast into the formulae Q4.19|) and (4.20). 



(Q Q Q^Q Q^) 
v ^234' ^134 ^124^123' ^1234 ' 



(4.21) 



these determi- 



In this way, we fulfilled the task of this subsection: all T-functions inside the T-hook 
can be expressed through 8 (ratios of) Q-functions. 

The whole procedure is illustrated in the attached Mathematica file, which for instance 
computes the expression of any Q-function in terms of the basis B\. 

24 In terms of the gauge fixing, decreasing the number of independent functions from 9 to 8 corresponds to 
adding the gauge constraint Q0 = 1. It leaves another gauge degree of freedom yet unfixed : corresponding 
to a gauge transformation Qj — > (J- nb -n f g) Qi, where rib and n/ are defined in footnote If we fix it, 
imposing for example Qi = 1, the number of independent functions becomes 7. 

25 Symbolically, we can still use ( |4.7| ) and (4.8), together with the rule 
Pi = -1, * e {1,2,3,4} and p Y = 1, IS {1,2,3,4}. 
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5. Asymptotic expressions for the generating functional and Q functions 



Now we will demonstrate our generating functional and Wronskian solutions of the 
AdS/CFT Y-system in the asymptotic, large size limit L — > oo. We will present the asymp- 
totic expressions of the relevant x~ an d Q-functions. All other Q-functions can be expressed 
through the basic ones through the Wronskian relations described above. Although it will 
be just a recasting of the ABA formulas of |6^] we feel that this Wronskian formulation is 
a right step in the direction of the derivation of the finite L system of FiNLIE's for the 
planar AdS/CFT spectrum. 

5.1 Asymptotic limit for the generating functional 

In the asymptotic limit the expressions for ^-functions entering the generating functional 



),(4.4) can be written explicitly in terms of the ABA Bethe roots. Namely, expliciting 



the Q-functions in the general expressions fl4.6|) we write pM |43| (see also p3[ ) 



B 3 R 1 B 3 R l Q2 B 1 R 3 Q2 B 1 i? 3 

B t R 7 c ^u B t R lQ&~ c ^XjQt +B 7 R 5 c ^u B 7 R 5 



where 



~F~ B^Rj ' B^RjQq ' Q(>B+R+ ' SJ B+R+ F+ 

(5.1) 

> ^ = ^T (5-2) 



and 



x 



H ={^) ^F S ^> H ={-) bT^W) (5 - 3) 

By A and A we denote the complex conjugations of A in the mirror and physical sheets, 
respectively (see the definition in |^|, §). 

Note that in the L — > 00 limit the first 4 £'s are suppressed in the mirror sheet whereas 
the last 4 are exponentially large. The expansion of the generating series can be interpreted 
as an expansion in the powers of a small exponential (which can be associated with the first 



factor in the first of the formulas (5.3)). It is then obvious that this expansion, which can 
be also qualified as an expansion in "wrappings" associated with the SYM planar diagrams 
wrapping around the operator, can be also understood as a series in powers of D = e~~' >: 



■2 



in the generating functional (|4.4|). The quantities R\ and are defined by 

R^\u) = l[(x{ U )- X l), B f\u) = ]l(^-xTj), KU7 (5.4) 

3=1 3=1 

where x t ■ = x{u h •), j ^ K\ encode the positions of Bethe roots. When the subscript / 
isn't specified, the value I = 4 is assumed, and by definition x { . = x{u l - ± |). The dressing 
factor is S(u) = FL <r(a;(ti), x 4 ■) where a is the BES dressing kernel [|^] (see |34], |6^] for a 
nice integral representation of the dressing kernel). 
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Furthermore, we define 

= B^R^ = | 11—^ | Ufa*- 1 < U 7 (5.5) 

5.2 Explicit expressions for the asymptotic Q-functions 

In the asymptotic (L — > oo) limit, where we know the expressions of all monomials £, 
all Q functions can be written in the leading ABA approximation in terms of Wronskian 
determinants. 

Matching the relevant ratios of Q-functions from ( f4.6| ) with the expressions fl5.1| ) we 
can express the Q functions of the basis B\ = {Qi, Q2, Q3, Q4, q\, q^^Qs^Qi} of the previous 
section in the asymptotic limit in the explicit analytic form given in the full generality in 
the appendix |B|. We demonstrate it bellow in a particular case of the SL(2) sector, which 
consists of the states such that V/ 7^ 4, Ki = 0, implying that Bi,R[ and Qi are equal to 1 
when / ^ 4. The basis of 8 independent Q-functions looks as follows: 

Q 1 = h(l-7) Q z = Q 1 .f- iu + iL±^\ (5.6) 

- 1 1 / + 

Qs = F ■-(!--) Q 4 = Q 3 - -*«+ (5-7) 



h F \ 2F-l r 

q i =[-M-^ F + ■ q % = f + hF (5.8) 

F Q + , Q~\ _ I 

ft 



where f,f,h and /i are defined by the recursion relations 

— o 

t f = H, ^ = H, £ = F, f — = F (5.10) 
n h+ 1 f 

while M is a polynomial function of u, defined by the recursion relation 

M+ - M~ = 2Q - Q (f~ + - Q ++ (f + + . (5.11) 



In the RHS of ( 5.11 ) we find a sum of two terms both having the form Q + ( F + = 



-q(+) [R^ + B( ^R^) but with different shifts of argument. It can be easily seen 

that they are polynomials 26 so that M itself is a polynomial which can be found (up to 



26 First, it is clear from ( ]5J ) that = Y[f=\ ^ is a constant. Then one can expand B (+) 7? (_) + 

B(-)flW = (n, x( u ) - 1^) (n, ^ - <J + (Uj x ( u ) - x Zj) (rij ^ - : the &st p^ct 

gives terms of the form (xj 4 , £4"^ , • ■ ■ , Jcjj ) (^J^ , ■ • ■ , X I j n ) x ( u ) n_ "\ while the second product gives 
terms of the form {x^f ,x^ i ,• • • , xj 4 ) (^4^1"' 1 j„) ^(u)" - *"- Then the definition of the Zhukovsky 
map- = x+l/x allows to rewrite each term (xj, , x^ i2 ,■■ • ,a:J im ) (^ji)'" i^Jjn) (z W I_m + x(tt) m ~ n ) 
as a polynomial in u. 
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an additive constant) from the eq,( 5.11| ) by matching the coefficients of the RHS and the 
LHS. 



For instance, for Konishi state, there are two roots (1(4,1, 1*4,2) = ( u i> ~ u i)i so that x; 



4,2 = "<i = Then 27 Q+ [F + ij = 2 x x 2 (4 - ^ + (x+ - ^ 



and ( pl| ) is solved by M = 2t«V ( - 5 ^~ + 8 + 2 ( x i ~ ^tJ ( x i 
which is indeed a polynomial in u. 

5.3 Physical symmetries 

Some symmetries can be identified in this asymptotic solution, which can be viewed as 
generalizations of symmetries of the characters of the classical monodromy matrix. 

Lower boundary: The first observation is that Qg = 1. This could already be immedi- 
ately seen from ( |5.1| ) by computing ^4- = (f^) (f-f-0 = ^' This is a generalization of 

the fact that the classical monodromy matrix £l(x) £ SU(2, 2|4) has the super-determinant 
equal to 1, and in terms of T functions, it means, by virtue of the last relation from (4.11), 



that To, s = 1 (as explained in Appendix |B.1| we can drop a possible i-periodic factor in 
%)• 

Mirror reality: The reality of the T functions on the mirror sheet, which is in the 



asymptotic limit a consequence of the relations [43] £1 = £2, £i = ^2' ^3 = £3 = £4' then 
translates into the condition that 28 




Q T = lit" 1 ) 1 fao-of ° 9i) Qe x (J) (5.12) 



which involves the gauge function g\ = -==, while nr and n& are the number of "bosonic" 

r t J 

and "fermionic" indices 29 in /, and the index exchange function / 1— > Ex(/), transforms all 
indices in / as follows: l«2,3f>4,if>2,3«4. 

For instance, when / = {3,4,2,3,4}, the equation ( 5.12| ) states that 30 



1 ^ 1 



Q 12i = Q34MI = + J x Q 43i43 = + - ^34134 (5-13) 

Wing exchange transformation Another symmetry of the asymptotic solution corre- 
sponds to the fact that, up to a gauge, T a ,_ s is obtained from T a , s by the transformation 
W£ : {-Bfc -B- -Bs-fc) Rk ^ R&-k}- In terms of Q functions, this relation reads 



V/, 



I "t> + "f I / \ 

Q / = (-l)»»(-l)L " J (J"„ b _ n/ o g 2 ) W£ (QwE(7j) (5-14) 



27 In these expressions, denotes x(u 4 1 ± ^). 

28 We use the definition (-if = (-1)\ 

29 n b = Card/ n {1, 2, 3, 4} and n f = Card/ n {i, 2, 3, 4} 

30 In (5.13), the second equality is given by (3.12), and the sign + comes from r[ig/(~l) 1 



( — )(+)(+)( — )(+) = +1. The third equality reorders the set of indices, and involves the sign of a per- 
mutation. This permutation exchanges the positions of two pairs of indices, hence the + sign in the third 
equality. 
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where the Wing Exchange function I i— > WE(7), transforms the individual indices in / 
according to 6 o 5 — b, /o5 — / = 5 — /, where 1 ^ 6 ^ 4 and 1 ^ / ^ 4 
and x h-> L X J denotes the floor function. For instance, for / = {1,2,1}, ( |5.14 ) states 



that 31 Q 12 j = — giYJE (^Q^jj = <?2W£ (Q12123)' which can be written as hf -B1-R3 



g2W£ [hf BjR^j, hence the value of the gauge function implied in that relation is 
92 = F 2 . 

5.4 Comments to the asymptotic Wronskian solution 

These equations simply recast (^J) hito the- form of determinants involving only the func- 
tions of the basis B\. The size of these determinants is fixed and doesn't increase with 
(a, s). 

On the other hand, the asymptotic Q-functions of B\ involve the functions / and 
h, which are expressed through infinite products and cannot be avoided in Q functions, 
because we have to express them knowing explicitly (in terms of £ functions) only the 
Q^/Qj . But fortunately in the T functions these infinite products are absent since the 

r 1 s | r_ s i 

T- functions are made of products of the type involving only finite products 

Similarly, although the Mjj functions look like infinite sums 32 it is possible to show 

\+s] f— si 

that T- functions involve only differences of the form M\- 1 — M\ ■ , which are finite sums. 
One can see it by certain subtractions of columns or lines of the determinants (4.1S- 4~15| ) 
which does not change their values. 

In conclusion, we demonstrated in this section our general solution of Hirota equation 
given in the previous section, in the case of asymptotic, large L limit of the AdS/CFT 
Y-system. As our experience with the principal chiral field model showed |jl9[ |20[| such 
asymptotic Wronskian expressions in terms of well chosen Q-functions can be very useful 
in establishing the finite L solution of the sigma model under investigation. 

6. Conclusion 

The main purpose of this paper was to express all the Y-functions and the associated 
T-functions entering the T-hook of the AdS/CFT Y-system, using its discrete integrable 
Hirota dynamics, in the form of an explicit Wronskian determinant expressions parameter- 
ized through a finite set of 8 Baxter-type Q-functions (7 independent Q-functions after all 
the gauge constraints are imposed). We view it as an important step towards the deriva- 
tion of a Destri-deVega-like finite system of non-linear integral equations (FiNLIE) for this 
important model. 

For certain relativistic sigma- models these Wronskian expressions, due to the integrable 
discrete Hirota dynamics of the Y-system, gave us a direct access to the corresponding 



The first equality is given by (5.14), and the sign — comes from ( — l) 2_l ~ L 2 J = —1. The second equality 



reorders the set of indices, and involves the sign of a permutation. 

32 In is not the case for instance in the sl(2) sector, because the M functions are polynomial and are easily 
identified, and no infinite sum arises. 
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FiNLIE |T|| ^Q j. In the case of SU(N) principle chiral field the corresponding Q-functions 
were simple polynomials in the asymptotic large L limit, and could be easily generalized 
to describe the finite size L system by introduction of certain N — 1 discontinuities, or 
"densities", vanishing at large L, along the whole real axis of the spectral parameter u. It 
suffices then to substitute these expressions into N — 1 TBA equations for the momentum- 
carrying nodes to write the needed FiNLIE |2C]. 



We don't know yet how to introduce these densities in the AdS/CFT case although 
it is conceivable that they could be non-zero only on the Zhukovsky-type cuts with the 
branch points at u = ±2g + y for some n € Z. Apart from the analyticity in the spectral 
parameter u, we also have to understand and incorporate into the Wronskian solution the 
asymptotic properties of T-functions w.r.t. large values of a and \s\. The hope is that, as 
in the case of SU(N) principal chiral field, there exists for the Wronskian solution of T- 
hook the "best" basis of 7 Q-functions with the simplest possible analytic properties. The 
analyticity of the rest of the quantities, Q-,T- or Y-functions, will be simply a consequence 
of the analytic properties of these 7 Q-functions and of the Wronskian formulas presented 
in this paper. We hope to describe it in the future work. 

Another interesting problem is to understand whether our Wronskian expressions could 
be promoted to an operatorial form. In the quantum spin chains or even in the confor- 
mal field theories the Q-operators enter the same commuting family of operators as the 
T-operators (transfer matrices), and operatorial form of Wronskian expressions makes a 



perfect sense and can be constructed, for example using the approach of [26, q6|, 57, 31] 



using (/-oscillator representations of the quantum affine algebra and the universal R-matrix. 

Let us also mention an interesting problem of the generalization of our Wronskian 
representation to another integrable duality, AdS4/CFT3, relating the 3D ABJM gauge 
theory and the sigma model on AdS^ x CP 3 . The Y-system for this model is a known 
|66|, [l4j solution and the corresponding Q-system was obtained for some cases in ]14| ]. 
The Y-system contains only one wing and it could be easier to study than the AdSs/CFT4 
duality. 

Apart from these main, physical tasks there are also a few other, more technical or 
mathematical questions to understand in our formalism. There should exist a natural gen- 
eralization of the Wronskian solution in the (2|4|2) T-hook given here, to all {M\\N\M2) T- 
hooks related to the infinite dimensional representations of u(M\, M2I-/V) The method 



of the "co-derivative" [67, 68] looks especially promising and could also be useful for op- 
erator construction mentioned above. This could potentially open an interesting field of 
research related to the search and investigation of a possible large new class of non-compact 
sigma models. 

The Wronskian determinant representation of the T-functions, solving Hirota equation 
in the AdS/CFT related T-hook in terms of a finite set of Q-functions, presented in this 
paper, gives us reasonable hopes for the construction of an AdS/CFT FiNLIE system, the 
analogue of Destri-deVega equations, and for a deeper understanding of the physical nature 
of AdS/CFT integrability. 
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A. Relations between Q functions 



These QQ-relations represent special versions of the general Pliicker relations for deter- 
minants. That means that there exists a certain number of QQ-relations expressing some 
Q-functions through the other, leaving only 8 independent functions. A few useful determi- 
nant relations express Q-functions of a later level of nesting through the ones on an earlier 
stages, such as 



Ql 



and 



(-1) 



m(m — 1) 



-m+n]\ 



^[-m+n+2fc-l]^ 
b 3 ) 



blb2...b m flf2---fn 



(Q[ _m +n] )V m _ ri _ 1 oQ, 



for m ^ n (A.l) 



b m flh---fn 



v n(n — 1) . / t \ 

(_]_) 2 + m ( n !) 


/ n [-m+nk 

/ n [-m+n-2i+l]\ 

' • Lfk )l^j^n—m,l^k^n 




m 

Fn-m-1 Ql 



for m ^ n. 

(A.2) 



For 7 = 0, these formulae reduce 33 to the determinant solutions of the QQ-relations in 
1 32]. Next, we introduce a useful trick on the index set for the determinant formulae, which 



33 In this sense, these formulae are a generalization of determinant formulae of Theorem 3.2 in [ p2[ . 
However, one can easily seen that these follow from Theorem 3.2 just by manipulating the index set 
of Q-functions. For two tuples, I and J = (6162 • • • 6,71/1/2 • ■ • /») (/ is fixed), let us consider a gauge 
transformation 

Qi,j 



Qi,j = 



(A.3) 



where '7, J' is a concatenation of / and J. Since there exists a relation Q 7 = Q 7 = 1, one can apply the 
Theorem 3.2 for Qj j. Then we obtain the formulae (A.l) and (A.2). 
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may be called "bosonization" or "fermionization" trick. Let us denote 

Qi = QlAB, Qi = QlAF 



(A.4) 



where AAC = (A U C) \ (A n C), B = (1,2,3,4) and F = (1,2,3,4). In other words, we 
define a Q/ through the corresponding Q j by adding to its indices forming the set J all the 
indices from the set B = (1,2,3,4) and then removing those of them which already were 
contained in J, to get finally I. 
From ( A.l ), we obtain 



0/ 



m(m — 1) 

-1)— s~ 



/A[-m+2fc-lk 



C\C2...C m 



F, 



m— 1 



(A.5) 



From ( A. 2 ), we obtain 



Qi 



-1) 



n(n — 1) 



,~{n-2j+l]s 



(A.6) 



This trick is efficient in the sense that it allows to write quite easily all Q functions in 
terms of only 8 functions, typically the single-indexed Q functions. 



B. Asymptotic expression of Q function 

In the general case (outside the sl{2) sector), the formulae (|5.6|)-(|5T9|) become, in full 
generality : 



Q 



Q 



3 



h {-FQ2 + Qt) 



B 3 R, 



Q 2 = Qi • I M 12 ' QlQ * FQ2 + Q * 



2Q2Q2 fq 2 ~ Qi 



> 2 v^ 2 i ^ 2 

F{-Q 6 + FQ { 

FhB 7 R 5 



Q 4 = Q 3 -(-M 34+ Q ^ Q « +FQ « 



ZQtQs - FQl 
g 2 = f+hFBrfz 

?4 = ~ 

ft 
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where M12, M34, Mj* and M-~r are meromorphic functions (without cuts, as explained in 
Appendix |B.1|) on the complex plane defined by the recursion relation 



12 12 




HTO Hit J J 

-Q++F+ + 2Q^ 

M+ - Ml* = — r — + , , Q2 (B.ll) 



2Qi Q3 



M+-ilC;=— ^ + — IU<^. (B.12) 

Note that all these recursion relations could be solved analytically in terms of infinite 
products, or of an integral representation taking into account their analyticity properties. 
In the sl(2) sector of quantum states of the AdS/CFT system, they can be identified simply 
as polynomials. 

We can check, by the direct substitutions of the formulas ( |B.1| )-( |B^ ) into the ex- 
pressions (4.6), that they reproduce (|5.1| ), and hence the right asymptotic solution of the 



Y-system given in 



B.l Hints of derivation of these asymptotic Q functions 



Let us explain the derivation of (B.1)-( B.8 ), starting from the expressions of the ^-factors 
(|5.l| ). We will see that the derivation assumes that a certain reality condition is satisfied. It 
is explained here only for the asymptotic limit (because ( |5.iD gives an explicit expression to 
start with), but in fact it applies to the construction of any real solution of the Y-system 34 , 
in terms of 7 independent Q-functions. 

The first step is to rewrite ( |5.1| ) in terms of Q function (by matching ( |5.1| ) with ( |4.6j )). 
This gives Qj = 2%^, Qjj = -Q h = Q 2 J, Q 12i = Q il2 = hf ~B\R 3 , Q 12i2 = Q il22 = 

/' Ql2123 = Qil223 = ^-^7^5/ j Ql23i23 = — Qil2233 = fQs> Ql234i23 = ^1122334 = 

-° ' dI1U W ~~ W 11223344 — L - 
h 

In particular it gives Qg = 1. This parameterization is done up to 8 arbitrary i- 
periodic functions Vj(u) = Vj(u + i), j £ I : the £ functions do not change under the 
transformation Qj (->■ Yiiel uj Card Q/ (where Card(I) denotes the length of I), which 
leaves the QQ-relations unchanged, and leads to an i-periodic gauge transformation of the 
T- functions. 

Then the QQ-relation Q is used to find Qi = (Q+ - Q~j)/Q£ = ^Rl (Qt ~ FQ2 )• 
The same QQ-relation also gives 



In the sense that (5.12) is satisfied for a gauge g\ which is a priori unknown. 
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Q 



Q+ 2i Qr-Q7 91 -Q^ h+h-F- 



12 



'121 



Q 



11 



Q>2 



B+R+' 



Then the "bosonic" QQ-relation ( f4.§| ) allows to write 
plugging here the previous two expressions we get 35 



Qt + ) 



+ 



^3 



■ Ch 
F~ 



(B.13) 
By 



QtQi 



Q2 

Qi 



02 

Qi 



Which is solved as 



A+ - A' 



A = — 



QiQz 



FiQzY-QzQt 



A 



Q1Q3 



Q2Q2 
(B.14) 



Q2 

Qi 



M12 + - ( A - A 



M+ 



M 



12 



-(A + + A + + A~ + A' 



(B.15) 



These two first steps give (B.l) and (B.2). The same manipulations help to find Qj 
and Qg : First, Q~- is extracted directly from the ^-functions, and then the "fermionic" 
QQ-relations give successively Q^gj and Q-^r. Then is extracted from the "bosonic" 

QQ-relation (O) 



Q0-Q- 

q+q- 1 which is solved in exactly the same manner 



as the equation for 3^. 



The fact that M-^r is meromorphic is nontrivial, and is implied by the polynomiality 
of A = Op- + FQ + = Q + (^F + entering there. This polynomiality is explained in the 
footnote and it essentially uses the fact that ^p- = = E.^ which is a consequence 

of the asymptotic level matching condition n^i = 1- This proves the meromorphicity 

x j 

of A % . 

Finally, the expressions of Q3,Q^,Qj and are obtained simply by the use of the 
complex conjugation transformation on Q functions. 
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